o 

(N 
o 

o 



cr. 

u . 
bJO, 



ITP-UU-ll/38, SPIN-11/29 

Symmetry breaking and the Goldstone theorem in de Sitter space 

Tomislav ProkopeC 
Institute for Theoretical Physics, Utrecht University, 
Postbus 80.195, 3508 TD Utrecht, The Netherlands 

We consider an 0{N) symmetric scalar field model in the mean field approximation 
and show that the symmetry can be broken (down to 0{N—1)) in de Sitter space. 
We find that the phase transition is of a first order, and that its strength depends 
non-analytically on the parameters of the model. We also show that the would-be 
Goldstone bosons acquire a mass, effectively becoming pseudo-Goldstone bosons. 
Our results imply that topological defects can form during infiation. 
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I. INTRODUCTION 

Ever since Kirzhnits and Linde [l| pointed out that thermal radiative effects can induce phase 
transitions in the early Universe, phase transitions have played a central role in the early Universe 
cosmology. In particular, they have been used to drive out-of-equilibrium phenomena which can 
lead to creation of the matter-antimatter asymmetry, preheating, formation of topological defects, 
etc. The effects induced by particle creation in an expanding Universe setting are quite delicate 
and have not yet been fully understood, albeit there is a large literature on the subject, a non- 



representative sample includes Refs. 
interacting theory. Ford and Vilenkin 



17|. Based upon a mean field analysis of a scalar self- 
19| pointed out a long time ago that the infrared effects 



in de Sitter space may restore symmetries spontaneously broken by the vacuum. However, a crucial 
sign error invalidated the conclusions of Ref. 18|]. The infrared effects in (quasi-)de Sitter spaces have 



received a considerable attention in recent literature, and several papers have been published 



20 



24| which have - just as Ford and Vilenkin - treated the problem in the mean field approximation 



fsee also Refs. 



25 



26| for some recent mean field results on fiat space). In particular, the above- 



mentioned sign error has been corrected in Refs. 



23 



24j | . The current consensus is that the infrared 
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effects in de Sitter space are strong enough to restore the broken 0{N) symmetry. While this is 
correct if the corresponding effective action is averaged over large (super-Hubble) distances, from 
the observational point of view the more relevant question is whether the symmetry gets broken 



or restored when averaged over some fixed physical scale [3l|- In this work we take the point of 
view that the effective action should be averaged over some fixed physical scale and we show that 
symmetries are then generally not restored in de Sitter space. We also give a simple criterion for 
symmetry restoration. For simplicity, we consider here only the global 0{N) symmetric scalar field 
model. 

For pedagogical reasons we begin in section [TT] by analysing a real scalar field (0(1) model). The 
central part of the paper is section UTTl where we analyse an 0{N) symmetric model on de Sitter 
space. Section [IV] is reserved for a discussion, and the Appendix for technical details on the de 
Sitter space propagator of a scalar field. 

II. A REAL SCALAR FIELD 



The free action of a real scalar field 0(x) in D space-time dimensions reads, 
S[(l>] = [ d^x^ 



where mo and A denote the field mass and quartic coupling, respectively, g^^ is the metric tensor, 
g'^'^ its inverse, and g = det[g^^]. The metric signature we use is (—,+,+,..). The vacuum is 
given by 0^ = for rrig > and 0^ = — Gmg/A when m-Q < 0. In the case when 0^ > the Z2 
symmetry (0 — )■ —(f)) of the action ([T]) is (completely) broken by the vacuum in which 0^ = —6171^/ A. 
If one goes through a quench from a high temperatures {e.g in an early Universe setting), when 
?T7,Q < domain walls form by the Kibble mechanism such that the state spontaneously breaks 
translation invariance. By causality, at least of the order of one domain wall forms per Hubble 
volume. Once formed, their energy density scales as oc 1/a^ (a denotes the scale factor), such that 
in decelerating spacetimes domain walls will dominate the energy density at late times (a menace 
to get rid of), while in accelerating spacetimes (such as inflation) they get diluted. 

Here we shall perform a mean field analysis, for which the effective potential (up to two loop 
order) is of the form, 

Vmf = Vo + ]^ml[(\)^ + «A(x; x)) + ^ (0^ + 60^A(x; x) + 3[zA(x; x)f) + ^Tr ln[zA(x; x)\ , (2) 

where rn^ denotes a bare mass term, A > a quartic coupling, 0(x) = (fi|0(x)|fi) is a mean field, 
is a state, and «A(x; a;') = (fi|T[50(a;')50(a;)] is the Feynman propagator for the field fluctuations 
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5<^{x) = — where T stands for time ordering. For simplicity, we have assumed that gravity 
is nondynamical and that all quantities (pix) and iA{x] x) are either constant or adiabatically varying 
in time (on a de Sitter background). Varying the mean field action 



MF 



.A] 



d^x\/ —g{x) 



A) 



with respect to and «A(?/; x) results in 



A 



— ran ^A('a;; x 



"^0 - ^(0^ + «A(x;x)) 



(3) 

(4) 
(5) 



where = {—g)^^^'^^^g^''^/^^JJ denotes the d'Alembertian (D^. = g^^V as it acts on a scalar 
quantity. Since the coincident propagator «A(x; x) is in general divergent (see the Appendix), these 
equations need to be renormalised. From the structure of equations (HHS]), a mass renormalisation 
suffices (32[. When viewed as a function of the ultraviolet cutoff A, the coincident propagator 



exhibits power law divergences oc A ,A , etc. (25|, |26|], which are automatically subtracted in 



dimensional regularisation. In the special cases when the spacetime dimension is even [D = 2n, 
n = 1,2, . . .), there is a logarithmic divergence in A, that manifests itself as a simple pole in the 
coincident propagator, iA{x;x)div oc 1/{D — 2n) (ra = 1,2, . . .), see Eqs. fl58l - [62|) in the Appendix. 
When this divergence is absorbed in the bare mass mg, one gets a finite, renormalised mass term. 



= mo + «A(x; a;)div, «A(x; x)div 



2 / 



D-1 



(6) 



47]-(D-l)/2 

where can be either positive or negative. Hence, the renormalised, manifestly finite, form of 
Eqs. (HHi is. 



-9 



Dx -w? - ^«A(x; a;)fin 



- ^(0^ + «A(x;x)fin) 



A 







(7) 



zA(a;;x') = i5 {x — x') 



where zA(x; x)fin = «A(x; x) —tA{x; x)div is the finite part of the coincident correlator, cf. Eqs. (jSH 
Since is (by assumption) slowly varying, we have □( 



m + -iA{x;x){in + -(p 
2 



0, and Eq. §^ yields 
/) = 0. 



This is solved by = or, when 0^ > 0, by 



— 3tA(x: x)f]r, > 



(9) 



(10) 
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At early times iA{x]x){i^ grows as given in flBOHB^ . reaching at late times the de Sitter invariant 
limit fl59|) . and the condition ([9]) becomes, 

6 = L^-J- > fill 

where m'^-p = d'^VuF / d(f)'^ is the mean field mass satisfying the mass gap equation 

X/ N f — 2m^ — A«A(x; x)fin, if 0^ > 



2 V ^1 m2 + |^A(x;a;)fin, if 0^ = 



Note that the mean field mass muF is also the mass of field fiuctuations in the correlator equation 
When Eq. (159|) is inserted into (fT2l) one gets, 



This is solved by. 



A/f^r (^) 



We see that - when the symmetry is broken, 0^ > - there is a minimum \vp?\ for which the gap 



equation ( !T4l) permits a meaningful (real) solution 33| : 



\vr?\>ml = \\ (15) 



In D = 2,3 and 4, m^^ = v/A/(47r)/7, ^f\j2H^I^ /-n, and v/3A/2i7V(27r), respectively. There is 
a simple way to determine the physically correct sign in Eq. (fT4|) . In the limit when H ^ 
the infrared enhanced fiuctuations are absent, such that one should recover the tree level mass, 
m^p — > —2rn?. This then implies that the physical mean field mass corresponds to the positive 
branch in ( 1T4|) . 



m^F = + - <^ (r > 0) . (16) 
On the other hand, when the symmetry is unbroken, 0^ = 0, Eq. ( lT2l) implies, 

mtiF - mV^p - ^m^r = (17) 

which is solved by. 



^MF=^ + V^ + ^ (r = 0). (18) 

Here we have dropped the solution with a negative sign in front of the square root, because for that 
solution m^F < 0; which is unacceptable on physical grounds (in this case the de Sitter invariant 
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state would be unstable under small perturbations). When the solution ( ITS]) is inserted into Eq. f[T^ . 
one can show that the only real solution for is = 0, consistent with the assumption of unbroken 



symmetry made in deriving Eq. (fT8|) 34 1. 



To summarise, we have found that, when 



XHDT (^) 

m < —m^^ ~ ~\l — 2 (-D+i)/2 — (vacuum breaks the Z2 symmetry) , (19) 



the infrared fluctuations on de Sitter space may not be able to restore the broken Z2 symmetry 
(0 — )■ —(f)) of the vacuum of a real scalar field ([T]). In this case 0^ = 3m^p/A > and Eq. ( |T6l) 
applies. Otherwise, when > — m^j, the Z2 symmetry is restored and Eq. f fTSj) applies. From 
Eq. (fT6|) we see that, in the broken symmetry case, there is a minimum mean field mass, given by 
(m^p)cr = — = m^j,, implying that for any finite coupling A, as |m^| increases from |m^| < m^j, 
to |m^| > mlj., the order parameter 0^ = 3m^p/A will experience a jump. 



^'^ - ^ - V 2Avr(^+i)/2 ' 
implying that the (phase) transition is of a first order. Notice that increasing H and decreasing A 
strengthens the transition. 

The above analysis has been performed under the assumption of adiabaticity, i.e. the time 
and spatial derivatives of have been neglected. Let us have a closer look at the validity of 
this approximation. When < — m^j., the Kibble mechanism necessitates formation of domain 
walls, which spontaneously break the de Sitter symmetry. The thickness of the domain walls is 
-^^dw ~ 1/"^MF5 implying that, when Ldw > ^/H (L^w < ^/H) the domain walls will be of a 
super-Hubble (sub-Hubble) thickness. In the former case the adiabaticity assumption for 0(x) and 
zA(x; x), which was used in the derivation of the mean field results above, may be unjustified because 
of the large spatial gradients and time derivatives of the fields that can in that case develop. Indeed, 
from the mean field equation (^^, which in general homogeneous expanding spacetimes reads, 

A 



+ 3Hdt 5- + m^F 



3V^' = 0, (21) 



we can estimate the spatial gradients of the fields as (V/a)0 ~ mMF0 (an analogous estimate holds 
for the time derivative, dt4> ~ itiuf4'), implying that one cannot anymore neglect the derivative 
terms in (l2Ti) . such that - in this case - a more sophisticated (numerical) analysis is required. 
Nevertheless, we can make a rough estimate of the critical mass as follows. When L^w > ^/H 
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the superluminal expansion of the space on super-Hubble scales will win over the forces which 
hold the domain walls together, such that they will stretch together with the expanding de Sitter 
space, restoring the symmetry. This simple argument suggests that the actual critical mass-squared 
(relevant in the full dynamical mean field setting fl2T]) ). above which the de Sitter symmetry gets 
restored and below which it gets broken, is given roughly by —mlj. ~ —H^, and not by Eq. ( fTOll . 
To get a more precise value however, a more complete (numerical) analysis of Eq. (^T^ is needed. 



III. THE 0{N) MODEL 



s 



d^xy/—g 



(22) 



We shall now consider the symmetry breaking in an 0{N) symmetric scalar field theory in the 
early Universe setting. Recall that this model allows for formation of (global) cosmic strings (when 
N = 2), global monopoles (when = 3), global cosmic textures (when = 4), etc. The free 
action of an 0{N) symmetric scalar field reads, 

N 1 ^ A ^ 

2^ J2'^dMi9u(i)a) - 2"^S '^^ ~ 4iv 



a=l a=l '- a=l -' 

Similarly as in the Brout-Englert-Higgs (BEH) mechanism, when mg < and A > the vacuum 
breaks the 0{N) symmetry to the 0{N — 1) symmetry, such that the resulting vacuum manifold 
is the N dimensional sphere, S*^ ~ 0{N)/0{N — 1) 27|. As a result, one of the scalars acquires 
a mass, while the other — 1 scalars remain massless. These massless excitations are known as 
Goldstone bosons. Unlike in the simple 0{N) model ( l22l) . at low energies the Goldstone bosons 
in the BEH mechanism acquire a mass and become the longitudinal excitations of the and Z 
bosons. For that reason they are known as pseudo- Goldstone bosons. We shall now see that the 
Goldstone bosons of the 0{N) model in de Sitter space (more generally in inflationary spacetimes) 
become massive due to the infrared (super-Hubble) enhancement of scalar correlations, and in that 
respect they can be considered as pseudo-Goldstone bosons. 

The mean field (two loop) effective potential of an 0{N) symmetric field (pa corresponding to 
the tree level action fl22|) reads. 



MF 



N 



a=l 



+ 



2^ 

AN 



N X 2 / N \ ^ 

a=l ^ ^ a=l ' h=\ 



bb{x;x) 



(23) 



N , N X 2 Af 

+ 4 ^ 0a0b«A„6(x; a;) + f ^ iAaa{x\ x) j + 2 ^ [it\ah{x\ x 

ah=\ ^ a=l ^ ah=\ 



+ -Trln(«A, 



{x\ x) 
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resulting in the following two loop effective action 

N 



s 



MF\<Pa, ^bc\ 



d^Xy/—g 



N 



(24) 



a=l a=l 

Varying the action with respect to (pa{x) and Af,a{x']x) gives the following (mean field) equations 
of motion (c/. Eqs. dMD), 



A 



N 



- ml - + tAbbix; x)) 



b=i 

N 



2A 



N 



(j)a{x) - — ^tAab{x;x)(t)b{x) = 



- ml - — ^ {(pl + 2Acc(x; x)) 



c=l 



tAabix; x') 



(25) 



b=i 



2A 



N 



[0a(a;)0c(a;) + «Aac(a;; x)] «Acb(a;; x') = Sab 



t6 {x — x') 



(26) 



c=l 



where 



^A,fe(x; x') = (n\T[8ct)b{x')5(t)a{x)]\nj (27) 

denotes the time-ordered scalar field propagator for scalar field fluctuations, 64)a{x) = 4)a{x) — 
{Vt\(j)a{x)\VL) . Just as in the one scalar case, Eqs. (l26l) can be renormalised by absorbing the infinite 
part of the coincident scalar propagator (!59|) into the bare mass mg, 

{N + 2)\H^'^T (^) 



2 2 I 

m = m^ + 



AT 



47]-(£'-l)/2 



-^(D-1)-^ 1 



(2^ 



resulting in manifestly finite renormalised equations analogous to Eqs. dZHS]). 

In the limit when the fields are slowly varying (D^a — 0), the renormalised form of Eq. ([25 
yields the following criterion for symmetry breaking, 



A 



N 



2A 



N 



m 



+ -^^(06+ lAbb{x; X)fin) (l)a{x) + ^ X] ?Aaf,(x; x)fin</)fe(x) = . 



N 



(29) 



The field mass matrix is obtained by taking a second field derivative of Vmf (or equivalently by 
taking a single derivative with respect to zA^a), 

N 



Mlb 



A 



m 



c=l 



2A 

^afe + ^ [0a0fe + ^Aab{x■, x)fin] 



(30) 



Notice that this result can be read off also from Eq. f l26|) . as is the mass term of the propagator 
lAab- Both Eq. ( 129|) and (!30|) contain in general off-diagonal terms. One can diagonalize them by 



an N X N dimensional orthonormal matrix R = {Rah)-, R ■ R^ = I, for which 





ab<Pb 



vV 



lAf-^ 
d 

^22 



lA^^ 



v 












^A^^ 



(31) 



Because of the unbroken 0{N — 1) symmetry, %Af^ are all equal for 2 < i < N. The diagonal form 
of Eqs. ([29H30]) is: 

A 



m 



+ - (0' + 3tAi,{x- a;)fin + {N - l)tAi^{x- x)fin) = . 



(32) 



for 2 < i < AT) and 



Mf = (Mfi 



d \2 



X 



m 



+ - [302 + 3«A?i(x; x)fin + (A^ - l)tAi^{x; x)fin] 



N 
X 



Ml = {M^y = m' + -[<f)' + zAi.ix; x)fi„ + {N + l)tAUx; x)fin] , {2 < t < N) . (33) 



Eq. 



implies that the 0{N) symmetry is broken when 



d\2 



r = (0?) 



A 



- 3tAl^{x; a;)fin - (A^ - l)iAUx; x){,^ > . 



(34) 



Otherwise, 0f = and the symmetry is unbroken. When this is inserted into fj33|) . we get that the 
mass terms (in the broken phase) become. 



^ (Mfj2 = ^02 = -2m' - ^ [3^Ai,ix■, x)^^ + (N - l)tAi,ix; x)fi,] 



d\2 



N 
2X 



[lAi (x; x)fi„ - 2 Af , (x; x)fin] , (2 < z < A^) 



(35) 



In the special case when A^ = 1 the first equation agrees with Eq. f|T2l) (provided, of course, one 
rescales the A in (!35l) as A — )■ A/6). 

With this, the renormalised and diagonalised form of Eq. ( l26l) becomes. 



x: X 



i6 {x — x') 



[a - M|] lA 



jj(x, X 



i6 {x—x') 



{2<i<N). (36) 



The implied stability of de Sitter space then demands that both {Mf^)'' = > {N > i > 2) and 
> 0. Next we insert the coincident propagator ( l59l) into ( l35l) to obtain: 
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Notice that positivity of implies that > > 0. The fact that the Goldstone bosons become 
massive on de Sitter background is reminiscent of the BEH mechanism, in which the Goldstones are 
'eaten up' by the longitudinal degrees of freedom of the and Z bosons, thus becoming massive. 

Equations f l37|) are the main result of this work. In order to analyse them, it is convenient to 
work with the following dimensionless quantities. 



(_2m2)' ^9 (_2m2)' iV7r(^+i)/2(-2m2)2 ' ^''''^ 



_ Ai^^r(^ 



after which Eqs. fl37j) become 

^ - (it " ii) • ''^'"i^-^y 

Before we perform a general analysis of these equations, notice that in the case when N = 1, the 
second equation decouples, and one gets 

fii-fj.j + 3XD = 0, (40) 

whose (physical) root (ssI] is 

^^l = l + ^\-3Xn. (41) 

The minimum critical mass is then determined by < (Ad)^ = 1/12 (1 > /uf > (/ii)cr = 
1/2), which accords with Eg. f|T5l) (when one takes account of the different definition of A in the 
Lagrangians and (122!) ) Sq . 

In the general case the gap equations fl39|) admit a small coupling expansion. Similarly as in the 
thermal case, the expansion parameter is y/XE, and hence non-perturbative, 

.2 _ . ... .^ ^ (iV-l)(4iV-21) ,/2 , ,,,,2 



K = 1 - (iV-1) VAz, - ^A^ - ^ '-^ '-X%' + 0{Xi) 

,2 _ ^ 4iV-5,3/2 , ,^,,2 



In figure[T]we show the rescaled masses ii\ and /i^ as a function of Xd- Both, the small coupling series 
solutions P2I) (dashed) as well as the full solutions of f l39l) (solid) are shown. We have plotted the 
heavy (Higgs) scalar mass /i^ = Ml/{—2m'^) (the upper curves starting at one when Ad = 0) and the 
pseudo- Goldstone mass /i^ = Mg/{—2rn?) (the lower curves starting at zero when Ad = 0) for = 1 
(the most extended green curves), N = 2 (the intermediate red curves), = 4 (the intermediate 
blue curves) and = 10 (the most squeezed gray curves). While the approximate solutions fH2|) 
continue for a while longer, the exact solutions end suddenly at a critical point, at which a minimum 
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■D 



0.00 0.02 0.04 0.06 0.08 

FIG. 1: The rescaled masses /if = Ml/[—2m?) (upper curves) and /x^ = Mg/[—2m?) (lower curves) as a 
function of the rescaled (dimensionless) coupling A/? defined in (j38|) . Both the exact solutions (solid) of 
Eqs. (|39p and the approximate solutions (j42|) (dashed) are shown. When viewed from the left to the right, 
we have plotted the cases: = 1 (green), N = 2 (red), = 4 (blue) and = 10 (gray). 

(critical) mass is reached. Just like for the case when = 1, where the mass parameter decreases 
monotonically from fif = I (at A/j = 0) to {fifjcr = 1/2 (at \d = {Xd)^ = 1/12), see Eq. ( BTj) . for 
N > 1, fil evolves monotonically from /if = 1 (at = 0) to some (/if)cr > 0, the end point being 
a function of A^. 

As we have seen in the analysis of the one scalar field case, this end point plays an important 
role as it tells us how the system behaves at the critical point (where the phase transition takes 
place). In order to study the critical behaviour in some detail, in figure [2] we plot the critical mass 
parameters (/if)cr (upper solid blue curve) and (/ipcr (lower dashed red curve) as a function of A^ 
(left panel) and as a function of {Xd)ct = max[AD] for a fixed A^ (right panel) (corresponding to the 
end points on figure [I]). One can show that, to good approximation, (A£))cr — S/[2{N + 2)^] (more 
precisely {Xd)ct — 1.62/[A^ + 2]^'^^), such that (A£i)cr approaches approximately quadratically zero 
as N ^ oo. A manifestation of this is the fractional power behaviour of /if ^ close to the origin for 
small (AD)cr (large A^) seen on the right panel in figure [2l 
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^ N 0.0 
20 0. 




FIG. 2: The rescaled critical (minimum) masses for the Higgs-hke excitation (/uf )cr = {Mf)cr/{—2m'^) (the 
upper sohd blue curve) and for the pseudo-Goldstones (/^^)cr = (-^^ )cr/(— 2m^) (the lower dashed red 
curve) as a function of (left panel) and as a function of the rescaled critical coupling {Xd)ct (defined as 
the maximum allowed Xd for a given N) (right panel). Individual points where is an integer are shown 
as gray dots. 

Just as in the one field case, the phase transition in the 0{N) model is of a first order. Indeed, 
from Eq. (152]) we see that 



N{M[ 



>0, 



(43) 



2A A 

with (/ii)cr plotted in figure |2l In fact, from the left and right panels on figure[2]one can read off that 
{nD,, ^ 2/[3v^] and oc {Xofd^ oc {\/N^f'^\m''\-^/\ implying that A^^ oc \m\N^ly\'^l\ 

which is to be compared with the single field result ( l20l) . where we found that (A0^)Ar=i oc A^^/^. 
Hence, in the large N limit the strength of the phase transition ( H3l) exhibits a qualitatively different 
dependence on A and m? than in the single field case f l20p . 

For completeness, we shall now briefly analyse the unbroken symmetry case. In this case (pf = 
(2 = 1, 2, .., A^) and Eq. ([32]) implies. 



+ A (3zAfi(a;; x)fi„ + {N - l)zA^2(a^; x)fin) > 



(44) 



and analogous steps as above yield. 
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Just as in the broken case, when = 1 we have, 



which agrees with Eq. (ITSl) . In fact, it is quite easy to obtain the general solution of equations (H5|) . 
Indeed, observe that the second equation can be written as, 

{^i\-^xl){\D + ^i\^il) = ^■ (47) 

The positivity of \di yUi and /x^ immediately implies that the only consistent solution is 

^ll = ^il. (48) 

It is not surprising that in this case all particles must have the same mass since the sjTiimetry is 
unbroken. With this, the first equation in (HSj) is easily solved. 



which can be also written as. 



,.2 ,.2 /^^ (A^ + 2)m2 

Ml = Ml = — + \j— + ^ ^ , (50) 

where rricr is given in f llSp . This generalizes the real field result f llSp to the 0{N) symmetric case. 
From Eq. (150|) one can easily see that Mf = Mg > 0, as it should be. Expanding the solution P9|) 
in powers of Xd, we see that - unlike in the broken case fj42l) - the result fHOl) is perturbative (and 
analytic) in A^, = = {N + 2)\d - 2{N + 2f\l + 0{\l). 

The same criticism of validity of the adiabatic approximation used in this section applies to the 
0{N) model as it did to the real scalar field model (see the end of section HB . 



IV. DISCUSSION 

We have analysed the 0{N) symmetric scalar field model fl22|) in the mean field approxima- 
tion (J23l) on de Sitter space. We have shown that symmetry breaking can occur, and that the 
would-be Goldstone bosons acquire a mass (see figure [I]) due to the enhanced infrared correla- 
tions in de Sitter space. Next we have studied the strength of the symmetry breaking transition, 
and we have shown that the (phase) transition is always of a first order. Curiously, the jump 
in the order parameter fH3|) exhibits a non-analytic dependence on the parameters of the model. 
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/S.cfP' oc |m|A^^/^A~^/'*, where \m\ and A denote the mass parameter and the quartic couphng of the 
model. 

While the mean field results are of their own interest, it would be desirable to investigate whether 
(and how) the mean field results presented here change when one includes higher loop corrections. 
A first step in this direction is taken in Ref. 2J| where the local contribution to the self-mass from 
the two loop (sun-set) diagram was estimated, and where it was found that, in the massless limit, 
the mean field mass-squared gets reduced by a factor 1 / \/2. 

Second, it is instructive to compare our results with the (old) stochastic theory results of Starobin- 
sky and Yokoyama j^], which is known to resum the leading log(a) corrections to infrared correlators 
on de Sitter space, see e.g. Q. From Eq. (23) of Ref. 28| we read (upon a rescaling, A — A/6), 

-Lch = ^(0^) = 1^^^^' - 0.1614V^i/^ , (51) 

which is to be compared with Eqs. ( ITSj) and ( |T5l) . which in the limit when m? — )■ and in D = 4 
yield m^p — )■ \/3Aif^/(47r). This then implies. 

Even though the difference in the results is modest, the question - which result is correct? - 
is, nevertheless, important. In the derivation of the stochastic result (15T|) . one assumes that the 
tree level potential remains unchanged, i. e. that for the late time behaviour the tree level potential 
should be used when stochastic theory is applied to infiation. At the moment there is no fundamental 
understanding concerning whether the tree level potential or some effective potential should be used 
in stochastic formalism. We close this discussion by noting that one can recover exactly the mean 
field result (m^p)m-^o from stochastic formalism, provided one replaces the tree level potential 
V = (A/4!)0^ by its Gaussian counterpart, V — > Vcauss = (-^Z^) ((0)^)0^- While this is suggestive, 
it does not ultimately tell us what is the correct answer. 
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Appendix: The de Sitter space propagator for massive and massless scalar fields 

Here we review some of the basic properties of the scalar propagator on de Sitter background in 
D space time dimensions. The time ordered (Feynman) propagator i/S. = obeys the equation, 

pD - m^] 2A++(a;; x') = i5^{x - x') , (53) 

where m is a mass and D^. = {—g)~^^'^di^i_g'^'^y/—gdu is the scalar d'Alembertian in D spacetime 
dimensions. The de Sitter invariance allows one to write the d'Alembertian in a de Sitter invariant 
form, 

m -y)-^,+D{2-y)--— zA{x; x') = -^=^6^{x - x') , (54) 



where y = a{ri)a{ri')H^[—{ri — t]')^ + \\x — x'W^] is related to the geodesic distance on de Sitter space 
i{x;x') as, y = 4 sin^(-?f£/2). Here a denotes the scale factor,?] is conformal time and x comoving 
coordinate. The unique solutions for the relevant propagators of the Schwinger-Keldysh (or in-in) 
formalism can be written in terms of the Gauss' hypergeometric function 2-^1 as follows, 

H^'' r(^ + z.^)r(^-.^)^^ ^-1 „.D., y-^ 



^A-^{x■x') = j^,^ ^- x.F.^^ + u^, — u^; I-"—), (55) 

where 

Here > represents the (renormalised) field mass parameter, which includes the renormalised 
mass, the mean field correction (the finite part of {X/2)tA{x; x)) and possibly also the term that 
originates from a nonminimal coupling, Am^ = ^D{D — l)if^ in the lagrangian, AC = —C,R4>^, 
where R = D{D — 1)H'^ is the Ricci scalar in de Sitter space. The functions y"^^ (a,/9 = ±) in 
Eq. ( !55l) denote, 

y++ = a{r])a{r]')H^[-{\r]-7]'\-teY + \\x-x'f] 
y^^ = a{ri)a{ri')H'^[—{ri — ri' + ieY+\\x — x'\\'^] 
y^^ = a(r])a{'r]')H'^[—{r] — 'r]' — ieY + \\x — x'\\'^] 

y— = a{r])a{r]')H'^[-{\ri-r]'\+ief + \\x-x'\\'^], (57) 
with e > infinitesimal. All propagators in ( 155|) have the same coincident limit. 
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Due to the last T function, this propagator exhibits a simple pole in even dimensions, D = 2, 4, 6, .., 
which reflects an ultraviolet (UV) logarithmic divergence. Of course, the leading UV divergence of 
the coincident propagator in de Sitter space is the same as that in Minkowski space, and it is of a 
degree D — 2, the subleading is of a degree D — 4, etc., the degree zero representing a logarithmic 
divergence. As it is well known, dimensional regular isat ion is blind to power law divergences (they 
are automatically subtracted by analytic extension), and exhibits only logarithmic divergences. 
The effect of the propagator fl58p can be considered in the weak curvature (Minkowski) limit, when 
^ (in which case one recovers the Minkowski space result plus small corrections) and in a 
strong curvature regime, in which ^ H^. Ignoring the UV divergence in (|58|) one can naively 
expand it in powers of rn? /H"^, and one obtains [37|, 



«A(a;; x) 



H^-'r (^) 

47r(D~l)/2 



D 



-ip{D-l)-^ 1 



D 
~2 



1 



-lE 



D-1 



+ 



9^(D+l)/2 ™2 ~^ \ M2 h 



(59) 



such that in D = 2,3,4 the 0{m-'^) terms are H"^ / {A-nrn^) , iJV(2vrW), and 3i^V(87rW), re- 
spectively. In our analysis in the main text we assume that both finite and infinite m-independent 
terms in (l59l) are absorbed in the physical definition of the mass term. 

The de Sitter invariant limit will be attained after some time during inflation. If the mass is 
very small (m^ ^ H'^), the propagator will at early times grow logarithmically with the scale 
factor (linearly with cosmological time). This can be seen by recalling that in the infrared ll|] the 
coincident propagator satisfies. 



(4vr)^/2 r (1 + vo) r (i - vo) 



D 

1 

2 



X 



\2Ha) 



Z)-l-2vo 



(60) 



where we took a(to) = = 1 and is an infrared (comoving) momentum cut-off. Notice that at 
early times (and in the limit when m — )■ and vn {D — l)/2) the coincident propagator grows 



logarithmically with time as (see also Refs. 29 



30|) 



zAfx; x) 



47]-(D~l)/2 

27r(^+i)/2 



V(^-i)-^ 1- 



D 



-lE + 



D-l 



(61) 
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This is to be compared with the Onemh-Woodard coincident propagator for a massless scalar field jo] 



[^^(-' -)]ow = ,^(._,)/. T{^) + 2.(^^W'2 ^ • (62) 

The logarithmic growth saturates when the propagator reaches the de Sitter invariant value flS^ . 
which characterizes the time scale at which the propagator (and thereby the state) becomes de 
Sitter invariant. 
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□ 

[31] As early as 1986 Ford and Vilenkin [19[] correctly reasoned: "[..] it would be completely consistent 

with all observations for our present Universe to be a de Sitter space with H^^ ~ 10^*^ yr. It would 
be very surprising if this were to have any effect upon symmetry breaking on terrestrial or subatomic 
scales." While this observation is correct, it does not follow from their analysis. 

[32] To renormalise the mean field effective potential ([2]), one also needs to renormalise Vq. Since Vq does 
not affect our analysis, we shall not renormalise Vq here. 

[33] A simple algebra shows that both mean field masses-squared m^p_|_ in (jl4p are positive, and moreover 
both are consistent with broken symmetry cp'^ > 0, see (jlip . which is to be contrasted to the conclusion 



in 



23|. 



23 



24i ]. where the sign error in [18] in front of A was 



[34] Equation ([T8]) has already been derived in Refs. 
corrected. 

[35] Recall that, in the limit when A^i — )• 0, the physical branch yields fif = M'^/{—2m?) = 1. 

[36] Formally, one can also solve Eq. (j39p for /x^ in the = 1 case, and one finds for the critical value, 

{fil)cr: = (\/T3 - 1)/12 ~ 0.217, which agrees with the results plotted in figures [1] and [2j 
[37] The 0{m?') term in Eq. (j59p has a divergent coefficient. To make the renormalised gap equation 

consistent, the divergent part of the 0{m?) term would have to be absorbed in the renormalised mass 

term m?. 

[38] Not surprisingly, the coincident propagator of a light scalar field ([60]) and the Onemli-Woodard coinci- 
dent massless scalar propagator (I6ip possess identical late time logarithmically growing terms oc ln(a). 
The time-independent parts do not agree, however. But this was to be expected, since these constant 
pieces do not have an independent physical meaning, as they can be absorbed in the mass counterterm 
of the self-interacting scalar theory. 



